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Abstract As the missing complement to existing studies of shell galaxies, we carried 
out a set of self-consistent N-body simulations of a minor merger forming a stellar 
shell system within a giant elliptical galaxy. We discuss the effect of a phenomenon 
possibly associated with the galaxy merger simulations — a presence of multiple 
generations of shells. 



Two- generation shell structure 

Galaxies with stellar shells are thought to be by-products of galaxy mergers 0. 
Most, though not all, previous models, e.g. l|5]|3]|2][4], relied on test-particle simu- 
lations. No systematic explorations of galactic models with these shell structures 
as merger debris via fully self-consistent N-body simulations, naturally involving 
the dynamical friction and the progressive decay of the accreted galaxy, were con- 
ducted. To bridge this gap, we decided to carry out a set of self-consistent simu- 
lations of a minor merger between a giant elliptical galaxy (gE with the mass of 
dark matter halo 8 • 10 12 M and stellar component 2- 10 11 M e ) and a satellite dwarf 
elliptical galaxy (2 1O 1O M in total), using the GADGET-2 code 0. 

In order to study differences in the resulting shell system formed in differently 
centrally concentrated mass distributions, we prepared simulations with the gE 
galaxy in two versions: a two-component Plummer and a two-component Hernquist 
model, with the same effective radius. The dwarf galaxy is then released on a radial 
orbit with initial velocity ~100km/s and distance of 200 kpc from the giant galaxy. 

In the first simulation, the core of the satellite passes through, returns and makes 
a second passage across the center of the primary galaxy (-~lGyr after the first 
passage). This event leads to creating the second generation of shells. To our know- 
ledge, this process has never been simulated in any previous study of the shell gala- 
xies, although predictions in this sense were made, e.g. JT). In the first approxima- 
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Fig. 1 The time evolution of the edges of three similar shell systems: results from the self- 
consistent simulation with gE galaxy modeled as the Plummer model (solid lines); from the test- 
particle simulation with the same initial conditions, but without the second passage of the satel- 
lite (boxes); and from the self-consistent simulation with gE galaxy represented by the Hernquist 
model (dotted lines). While the stars, gradually released from the potential of the satellite, oscillate 
in the potential of the gE galaxy, they are forming shells alternately interleaved on both sides of the 
gE galaxy along the merger axis — the first shell on the right side (R) is later followed by the second 
shell on the left (L), etc. The velocity of the shell-edge expansion — slope of d n (t) dependency — 
is given by the gravitational potential and by the shell ordinal number within a given generation. 
Therefore the first shells from each generation move with the same velocity. Such a shell galaxy 
arisen in the first self-consistent simulation, if observed in the particular time (e.g. 1.7 Gyr after 
the first passage), may appear to cause "the problem of a missing shell", since the third (originally 
fourth) outermost shell is then detected on the same side (L) as the second one. 



tion, we can look at this as a new collision between the returning core part of the 
satellite and the gE galaxy. Within the same generation, the shells of the debris sys- 
tem are moving with decreasing velocity. As the subsequent passage is not present 
in the latter simulation (with a two-component Hernquist model for the gE galaxy), 
the subsequent shells created after ~1 Gyr move with different velocities compared 
to those belonging to the next generation in the former simulation, see Fig.Q] 
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